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Abstract 



' A parameter free investigation of the moments of inertia of ground state 

: 

, rotational bands in well deformed rare-earth nuclei is carried out us- 

O ■ 

ing Cranked Relativistic Hartree-Bogoliubov (CRHB) and non-relativistic 
I Cranked Hartree-Fock-Bogoliubov (CHFB) theories. In CRHB theory, the 

relativistic fields are determined by the non- linear Lagrangian with the NLl 
force and the pairing interaction by the central part of finite range Gogny DIS 
O , force. In CHFB theory, the properties in particle-hole and particle-particle 

Ch ■ channels are defined solely by Gogny DIS forces. Using an approximate parti- 



cle number projection before variation by means of the Lipkin Nogami method 
improves the agreement with the experimental data, especially in CRHB the- 
ory. The effect of the particle number projection on the moments of inertia 
and pairing energies is larger in relativistic than in non-relativistic theory. 

PACS numbers: 21.60.-n, 21.60.Jz, 27.70. +q 

Keywords: Relativistic and non-relativistic mean field theories, approxi- 
mate particle number projection, moments of inertia, pairing energies, charge 
quadrupole moments. 

Typeset using REVT^ 



*on leave of absence from the Laboratory of Radiation Physics, Institute of Solid State Physics, 
University of Latvia, LV 2169 Salaspils, Miera srt. 31, Latvia 



1 



I. INTRODUCTION 



One of the oldest problems in our understanding of collective motion of nuclei is the 
moments of inertia of ground state rotational bands in well deformed nuclei. They depend 
in a very sensitive way on the collective properties such as deformations and on pairing 
correlations of these many-body systems. Since rotational bands have been detected in nuclei 
nearly fifty years ago and since the first microscopic calculations of the moments of inertia 
by Inglis these quantities have been used as a testing ground for nearly all microscopic 
theories of collective motion. They describe the response of the strongly interacting nuclear 
many-body system to an external Coriolis field breaking time reversal symmetry. They 
are, therefore, in some sense comparable to the static magnetic susceptibility in condensed 
matter physics. 

The earliest microscopic calculations were based on a mean field of a deformed harmonic 
oscillator In these calculations, residual interactions were neglected. In this way one 

founds the values of the moments of inertia identical to those of a rigid body with the 
same shape, in strong disagreement with the experimentally observed values, which where 
considerably smaller. It has been pointed out already very early 0,^] that residual two-body 
interactions would lower the values of the moment of inertia obtained in the Inglis model. 
The most important correlations causing such a reduction are pairing correlations In fact, 
Belyaev |PJ^ showed that a simple extension of the Inglis formula in the framework of the 
BCS theory is able to reduce the theoretical moments of inertia dramatically because of the 
large energy gap in the spectrum of quasiparticle excitations occurring in the denominator 
of the Belyaev formula. Therefore, the small moments of inertia of the rotational bands 
provided one of the most important experimental hints for a superfluid behaviour of these 
heavy open shell nuclei. Extended calculations using the theory of Belyaev have been carried 
out by Nilsson and Prior using the BCS model based on the single particle spectrum of 
the Nilsson potential. 

Apart from the fact that the results of these calculations were relatively successful, there 
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are, as we know today, a number of open problems, namely: 

i) Belyaev's formula is based on generalized mean field theory violating essential sym- 
metries. It was pointed out already by Migdal [§,0], that Galileian invariance is broken. 
He therefore modified the Belyaev formula by taking into account more complicated cor- 
relations to correct the violation of this symmetry. The question of the restoration of the 
broken Galileian invariance in the particle-hole and particle-particle channels has been later 
discussed in a number of articles, see for example Refs. |1TlHT3|] and references therein. 

ii) Since Belyaev's formula describes only quasiparticles moving independently, higher 
order correlations have to be taken into account. This has been done by Thouless and 



Valatin who considered all orders of the interaction in a theory describing the linear 



response of the system to the external Coriolis field. Marshalek and Weneser showed 
that the method of Thouless and Valatin preserves all the symmetries violated in the mean 
field approximation in linear order. In that sense Migdal's formula was just a special case 
to deal with Galilean invariance. Marshalek showed in a series of papers (see for example 
Ref. |T^) that this is just the linear approximation of a more general theory based on Boson 
expansion techniques treating the symmetries appropriately in all orders |16| ]. 

iii) Much more elaborated versions of the cranked Nilsson model [0,0 showed that the 
Z^-term in this model, which corrects in an elegant way the fact that realistic potentials 
for heavy nuclei are much flatter than an oscillator in the nuclear interior, introduces a 
strong spurious momentum dependence. This leads to the values for the moments of inertia 
deviating considerably from the experimental values. However, this problem is to large 
extent cured either by Strutinsky renormalization of the moments of inertia [1^ or by an 



additional term to the cranked Nilsson potential that restores the local Galilean invariance 
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Realistic applications of the Thouless- Valatin theory are by no means trivial. They 
should be based on self-consistent solutions of the mean field equations, because only for 
those solutions the RPA theory preserves the symmetries [PTI. In addition, they require the 



inversion of the RPA-matrix. Meyer-ter-Vehn et al. have carried out such calculations 



in a restricted configuration space replacing the self-consistent mean field in an approximate 
way by the Woods-Saxon potential. As residual interaction they used density dependent 
Migdal forces in the p/i-channel and in the pp-channel. These interactions have been 
carefully adjusted to experimental data for the underlying configuration space. The results 
of these calculations showed, that there are indeed effects originating from both channels, 
each of them modifying the Belyaev values, but cancelling themselves to a large extent. 
Therefore one could understand why older calculations based on the generalized mean 
field model gave reasonable results as compared to the experiment. 

Nowadays there are theories available where the Hartree-(Fock)-Bogoliubov equations 
can be solved in a fully self-consistent way in the rotating frame for finite angular velocity 
Q. Using the resulting wavefunctions l^^^) the Thouless-Valatin moment of inertia can be 



found as 

(1) 



J = 4?^{'^n\Jx\'^n) 
ail 



Q=0 

In this way one avoids the inversion of the full RPA-matrix, as task which is so far tech- 
nically impossible for realistic forces in a full configuration space. Among these theo- 
ries the properties of rotating nuclei are described in a way free from adjustable param- 
eters only in the Cranked Relativistic Hartree-Bogoliubov (CRHB) theory |p3ip4| and non- 
relativistic density-dependent Cranked Hartree-Fock-Bogoliubov (CHFB) theory with finite 
range Gogny forces p5| , p6| . Several realistic investigations of the moments of inertia in 
normal- and in super-deformed bands have been carried out in the literature in the frame- 
work of non- relativistic CHFB theory with Gogny forces p5| , p7| - p0|] . Similar investigations 



in the relativistic framework have been performed without pairing in the A ~ 60 [31], 80 



and 140 — 150 p3|-p5|] regions of super deformation where the pairing correlations are 



expected to be considerably quenched at high spin. The recently developed formalism of the 
CRHB theory has been applied so far only for the description of the moments of inertia in 
the A ~ 190 mass region of superdeformation [^,^. A very successful description of the 
moments of inertia has been obtained in the framework of these two theories. The aim of 



the present investigation is to find the similarities and differences between these two theories 
using in a systematic way the moments of inertia of rare-earth nuclei as a testing ground. 



II. THEORETICAL TOOLS 



The CRHB theory p3|,^ is an extension of cranked relativistic mean field theory 



36| , |33|j3^ to the description of pairing correlations in rotating nuclei. It describes the nu- 



cleus as a system of Dirac nucleons which interact in a relativistic covariant manner through 
the exchange of virtual mesons |^: the isoscalar scalar a meson, the isoscalar vector uj 
meson, and the isovector vector p meson. The phonon field {A) accounts for the electromag- 
netic interaction. The CRHB equations for the fermions in the rotating frame are given in 
one-dimensional cranking approximation by 

^][ ]=eJ ] (2) 
where h = hu — X is the Dirac Hamiltonian ho for the nucleon with mass m 

hD = aHV -V{r)) + Voir) + p{m + S{r)) (3) 
minus the chemical potential A defined from the average particle number constraint 

{^n\N\<Pn) = N. (4) 

The Dirac Hamiltonian contains a repulsive vector potential Vo(r), an attractive scalar 
potential S{r) and the magnetic potential V{r) which leads to non- vanishing currents in 



the systems with broken time- reversal symmetries These currents play an extremely 

important role in the description of the moments of inertia [^,0 and thus they are taken 
into account fully self-consistently in the calculations. In Eq. (|^), Uk and Vk are quasiparticle 
Dirac spinors and denotes the quasiparticle energies. Furthermore, and are the 
projection of total angular momentum on the rotation axis and the rotational frequency. 

The time-independent inhomogeneous Klein-Gordon equations for the mesonic fields are 
given by 



{-A - {QM' + ^1} ^(r) = -9. [p'sir) + p:{r)] - g2<y\r) - g,<y\r), 
{-A - {n^L^f + ml) oooir) = [ff,{r) + p^ir)] , 
{-A - {f2,{L, + S^yf + ml} a;(r) = g^ [fir) + j"(r)] , (5) 

with source terms involving the various nucleonic densities and currents 

k>0 k>0 

fir) = T.iyk{r)]^^nir). (6) 

fc>0 

The sums over k > run over all quasiparticle states corresponding to positive energy single- 
particle states [no-sea approximation) and the indexes i could be either n (neutrons) or p 
(protons). For simplicity, the equations for the p meson and the Coulomb fields are omitted 
in Eqs. (H) since they have the structure similar to the equations for the u meson, see Refs. 
|23| , |2^ for details. Since the coupling constant of the electromagnetic interaction is small 
compared with the coupling constants of the meson fields, the Coriolis term for the Coulomb 
potential Ao(^) the spatial components of the vector potential A{r) are neglected in 
the calculations. 

The pairing potential A in Eq. (^ is given by 

A^Aa. = iT^VaL^c, (7) 

cd 

where the indices a,b, . . . contain the space coordinates r as well as the Dirac and isospin 
indices s and t. It contains the pairing tensor k 

K = Kir,s,t,r',s\t) = J2v:ir,s,t)Ukir',s\t) (8) 

Ek>0 

and the matrix elements V^^d effective interaction in the pp-channel. In the present 

version of CRHB theory, pairing correlations are only considered between the baryons, be- 
cause pairing is a genuine non-relativistic effect, which plays a role only in the vicinity of 
the Fermi surface. The central part of the Gogny interaction containing two different finite 
range terms (see Eq. (^) is employed in the pp (pairing) channel. 
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The CRHB calculations have been performed with the NLl parametrization of the 



RMF Lagrangian. The DIS set of parameters is used for the Gogny force in the pairing 
channel. The CRHB-equations are solved in the basis of an anisotropic three-dimensional 
harmonic oscillator in Cartesian coordinates. A basis deformation of jSo = 0.3 has been 
used. All fermionic and bosonic states belonging to the shells up to Np = 13 and Nb = 16 
are taken into account in the diagonalisation and the matrix inversion, respectively. This 
truncation scheme provides reasonable numerical accuracy for the physical observables which 
as estimated in the calculations with larger fermionic basis is on the level of ~ 1.5% or 
better for kinematic moment of inertia J^^-* and charge quadrupole moments Qq. In order 
to calculate the derivative with respect to Q in Eq. (|l|), all CRHB calculations have been 
performed at rotational frequency = 0.05 MeV. 

The starting point of the non-relativistic CHFB theory based on the Gogny force is the 



phenomenological finite range two-body interaction of the form 

1=1,2 

+ iWLsiy 12 A 5(ri - r2)Vi2)(o-i + cts) 

+ h{l + P'^xo)5(ri - r^MR)^'' (9) 

which is used simultaneously both in pp and ph channels. In Eq. @, R = {vi -\- r2)/2. 
The transformation to the rotating frame leads to equations similar to (HJ^J^). The only 
difference is that the Dirac Hamiltonian of Eq. (|^) is replaced by the non-relativistic Hartree- 
Fock Hamiltonian hij containing the density-dependent Gogny force and the rearrangement 
term dTij, stemming from the density dependence of the force 

(6H' \ 
— (10) 

In the above expression fij{r) is the quantity appearing in the second quantization form of 
the density operator p(r) = fij{.f)ctcj. The parameter set DIS has been used in the 
present calculations. The CHFB-equations are again solved in the basis of an anisotropic 
three-dimensional harmonic oscillator in Cartesian coordinates with the oscillator length 



bo = 1.98 fm and the deformation of basis Po = 0.3. Only single-particle states satisfying 
the condition 

hUxflx + hUyUy + hUzUz < NmaxflUJo with N^ax = H-l (H) 

have been included in the basis. The HFB equation has been solved with the Conjugated 
Gradient Method PO . 



We also consider in this investigation the fluctuations in the pairing field by using the 
technique of an approximate particle number projection before the variation introduced by 
Lipkin and Nogami (further APNP(LN)) and discussed in detail in non-relativistic case in 
Refs. i2|,|l],|0|. In the relativistic case, the same approximate particle number projection 



is used but only the pp-pait of the interaction is taken into account for the Lipkin-Nogami 



procedure, see Ref. ||2J] for details. 



III. RESULTS AND DISCUSSION 

In order to see the dependence of the results on proton and neutron numbers several 
nuclei in the Gd, Dy, Er and Yb isotope chains, the ground state rotational bands of which 
are close to rotational limit (i?(4+)/£'(2+) ^ 3.3), have been selected for the present study. 
The results of relativistic and non-relativistic calculations with and without APNP(LN) are 
presented in Tables 1-3 and Figs. 1-2 and compared with the experiment. Such quantities 
as charge quadrupole moments [deformations] , pairing energies and moments of inertia are 
discussed below. 

The calculated and experimental charge quadrupole moments Q and quadrupole defor- 
mation parameters /? derived from Q by 



/Ifivr 3 

Q = J ——ZRlp where i?o = 1.2A^/=^ (12) 
V 5 47r 

are shown in Table 1 and Fig. 1. The general feature is that the charge quadrupole moments 
Q obtained in relativistic calculations are larger than the ones of the non-relativistic calcu- 
lations. In the non-relativistic case, the Q values calculated with APNP(LN) are slightly 
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smaller than the ones obtained without APNP(LN) because of the larger pairing correlations 
(see Table 2), which in general favours more spherical configurations. This trend also per- 
sists in the relativistic case, but there are the cases (i64,i66j^j^^ i64,i66y]3) in which APNP(LN) 
leads to larger charge quadrupole moments as compared with unprojected calculations. As 
shown in Fig. 1, the non-relativistic results are somewhat closer to the experiment than 
the relativistic ones. However, considering that the experimental values of Q are subject of 



considerable experimental errors [^, one can conclude that both theories describe experi- 
mental charge quadrupole moments reasonable well which allows us to proceed further with 
the study of more sensitive quantities such as moments of inertia. 

In Hartree-(Fock)-Bogoliubov calculations the size of the pairing correlations is usually 
measured in terms of the pairing energy defined as 

Epair = -^Tr(AK). (13) 

This is not an experimentally accessible quantity, but it is a measure for the size of the pairing 
correlations in the theoretical calculations. These quantities are shown in Table 2 and Fig. 
2 for protons and neutrons separately. Both in relativistic and non-relativistic calculations, 
we observe that APNP(LN) leads to an increase of the pairing energies. This increase shows 
large variations as a function of the proton and neutron numbers. In general, this increase is 
larger in relativistic calculations. For example, proton pairing energies increase on average 
by factor 2.16 (with minimal and maximal increase being equal to 1.53 (^^^Dy) and 3.36 
(^^'^Gd)). In non-relativistic calculations, the average increase of proton pairing energies is 
only 1.55 (with minimal and maximal increase being equal to 1.11 (^^^Dy) and 1.97 (^™Yh)). 
Neutron pairing energies behave in a similar way but there the difference between relativistic 
and non-relativistic calculations is smaller: the average increase of neutron pairing energies 
due to APNP(LN) is 1.81 in the relativistic and 1.73 in the non-relativistic calculations. 
In some cases, as for example in the Gd isotopes, the increase of neutron pairing energies 
due to APNP(LN) is larger in non-relativistic calculations. This increase of pairing energies 
due to APNP(LN) will lead to an increase of the pairing gaps, as it is well known from 
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many phenomenological calculations using the monopole pairing force We also see 

that with few exceptions the pairing energies are smaller in relativistic calculations. An 
additional effect of APNP(LN) is the increase of absolute values of binding energies. In 
relativistic calculations, APNP(LN) provides an additional binding by ^ —2.5 MeV. 

Calculated moments of inertia are given in Table 3 and Fig. 1. Comparing the results of 
calculations without APNP(LN), it is clear that the moments of inertia are systematically 
larger in the relativistic case. Although one cannot completely exclude that this feature 
is to some extent connected with a different angular momentum content of single-particle 
orbitals in relativistic and non-relativistic calculations, a detailed analysis of pairing energies 
and moments of inertia suggests that this fact can be explained in a more realistic way by 
the different effective masses of the two theories: m*/m ~ 0.6 in RMF theory and ~ 0.7 in 
the non-relativistic theory. Thus the corresponding level density in the vicinity of the Fermi 
level is smaller in the relativistic theory which in general leads to weaker pairing correlations 
(see Table 2 and discussion above) as compared with non-relativistic calculations and as a 
result to larger moments of inertia. APNP(LN) restores to large extent the correct size 
of pairing correlations and thus its effect is larger in relativistic calculations. The average 
decrease of the moments of inertia due to APNP(LN) over the considered set of nuclei is 
1.35 and 1.15 in relativistic and non-relativistic calculations, respectively. It is also clearly 
seen that APNP(LN) improves in average and especially in relativistic case the agreement 
between experimental and calculated moments of inertia. The level of agreement between 
calculations with APNP(LN) and experiment is similar in both theories, however, some 
discrepancies still remain. 

IV. CONCLUSIONS 

In conclusion, the moments of inertia, charge quadrupole moments and pairing energies 
of well-deformed nuclei in the rare-earth region have been investigated within relativistic and 
non-relativistic mean field theories with and without approximate particle number projection 
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by means of the Lipkin-Nogami method. With no adjustable parameters it was possible 
to obtain good description of experimental data. It was found that the particle number 
projection plays a more important role in the relativistic calculations most likely reflecting 
the lower effective mass. In addition, it has larger impact on the moments of inertia and the 
pairing energies as compared with the charge quadrupole moments. Remaining deviations 
from experimental data could be related either to the parametrization of the mean field or 
to the interaction in the pairing channel or to the approximate character of the particle 
number projection. Further and more systematic investigations are needed for clarification 
of the main source of discrepancies between theory and experiment. 

A.V.A. acknowledges support from the Alexander von Humboldt Foundation. This work 
has been supported in part by DGICyT, Spain under project PB97-0023 and the Bundesmin- 
isterium fiir Bildung und Forschung under the project 06 TM 875. P.R. wishes to express his 
gratitude to the Spanish Ministry for Education for support of his work at the Universidad 
Autonoma de Madrid in the framework of the A. von Humboldt - J.C. Mutis program for 
the Scientific Cooperation between Spain and Germany. 

V. FIGURE CAPTION 

Fig. 1. Experimental and calculated charge quadrupole moments Q (top panels) and 
moments of inertia 2 J^^^ (bottom panels) of well deformed rare-earth nuclei. The experimen- 
tal data arc shown by solid unlinked circles. The results of calculations with APNP(LN) are 
shown by the lines without symbols. The lines with open symbols are used to indicate the 
results of calculations without APNP(LN). Solid and dashed hues are used for relativistic 
and non-relativistic results, respectively. 

Fig. 2. Calculated neutron (top panels) and proton (bottom panels) pairing energies. 
The results of calculations with (without) APNP(LN) are shown by the lines without (with) 
symbols. Solid and dashed lines are used for relativistic and non-relativistic results, respec- 
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tively. 
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TABLES 

TABLE L The calculated and experimental charge quadrupole moments Q and quadrupole 
deformation parameters (3 [shown in squared brackets] for typical well deformed nuclei in the 
rare earth region. The results of relativistic calculations are indicated by 'CRHB', while the 
results of non-relativistic calculations by 'Gogny'. The calculations without and with approximate 
particle number projection by means of the Lipkin-Nogami method are shown in columns marked 
by 'without projection' and 'with projection', respectively. The experimental data are taken from 
Ref. |3|. 



without projection with projection Expt. 





A 
A 




CRHB 


Go 


gny 


CRHB 


Go 


gny 






154 


6 


715 


[0.335] 


n r\'i A 

6.074 


[0.303] 


5.907 


[0.295] 


5.606 


[0.280] 


6.221 [0.310] 




156 


7 


199 


[0.356] 


6.792 


[0.336] 


6.886 


[0.341] 


6.601 


[0.327] 


6.830 [0.338] 




158 


7 


383 


[0.362] 


7.077 


[0.347] 


7.262 


[0.356] 


6.961 


[0.341] 


7.104 [0.348] 




IbU 


7 


577 


[0.369] 


7.286 


[0.354] 


7.490 


[0.364] 


7.200 


[0.350] 


7.265 [0.353] 


Dy 


15o 


5 


860 


[0.281] 


5.994 


[0.287] 


5.610 


[0.269] 


cr /loo 

5.438 


[O.iolJ 


6.107 [0.293] 




158 


7 


032 


[0.334] 


6.990 


[0.333] 


6.711 


[0.319] 


6.538 


[0.311] 


6.844 [0.326] 




160 


7 


496 


[0.354] 


7.297 


[0.344] 


7.373 


[0.348] 


7.102 


[0.335] 


7.13 [0.337] 




162 


7 


711 


[0.361] 


7.492 


[0.350] 


7.697 


[0.360] 


7.382 


[0.345] 


7.28 [0.341] 




164 


7 


928 


[0.368] 


7.626 


[0.354] 


7.883 


[0.366] 


7.543 


[0.350] 


7.503 [0.348] 


Er 


164 


7 


671 


[0.345] 


7.585 


[0.341] 


7.791 


[0.351] 


7.522 


[0.339] 


7.402 [0.333] 




166 


8 


047 


[0.359] 


7.781 


[0.347] 


8.075 


[0.361] 


7.728 


[0.345] 


7.656 [0.342] 




168 


8 


213 


[0.364] 


7.838 


[0.347] 


8.151 


[0.361] 


7.831 


[0.347] 


7.63 [0.338] 




170 


8 


137 


[0.358] 


7.899 


[0.347] 


8.075 


[0.355] 


7.782 


[0.342] 


7.65 [0.336] 


Yb 


164 


6 


552 


[0.287] 


6.900 


[0.302] 


6.602 


[0.289] 


6.828 


[0.299] 


6.60 [0.289] 




166 


7 


339 


[0.318] 


7.594 


[0.330] 


7.653 


[0.332] 


7.508 


[0.326] 


7.19 [0.312] 




168 


8 


362 


[0.360] 


7.997 


[0.344] 


8.222 


[0.354] 


7.864 


[0.339] 


7.59 [0.327] 




170 


8 


546 


[0.365] 


8.147 


[0.348] 


8.354 


[0.357] 


8.013 


[0.342] 


7.57 [0.324] 
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TABLE II. Pairing energies for typical well deformed nuclei in the rare earth region. For 
details of this table see caption of Table 1. 



'pair pair 



without projection with projection without projection with projection 





A 


CRIIB 




CRIIB 




CRIIB 


Q( )gii\' 


CRIIB 


(-rogny 


Gd 


154 


-6.790 


-7.413 


-11.264 


-14.566 


-4.556 


-8.857 


-10.176 


-10.828 




156 


-6.039 


-6.236 


-11.552 


-12.871 


-3.635 


-7.740 


-8.977 


-10.778 




158 


-7.071 


-6.511 


-11.495 


-12.355 


-3.151 


-7.247 


-8.748 


-10.487 




160 


-7.174 


-6.731 


-11.385 


-11.961 


-2.579 


-6.830 


-8.686 


-10.341 


Dv 


156 


-5.965 


-8.017 


-10.904 


-15.610 


-7.293 


-9.849 


-11.167 


-10.888 




158 


-7.007 


-8.085 


-11.742 


-13.870 


-5.596 


-7.773 


-9.888 


-11.234 




160 


-7.706 


-8.067 


-11.712 


-13.101 


-3.077 


-6.980 


-9.470 


-11.090 




162 


-6.331 


-7.887 


-11.440 


-12.581 


-2.859 


-6.323 


-9.289 


-10.837 




164 


-5.237 


-6.755 


-11.046 


-12.169 


-3.497 


-5.952 


-9.187 


-10.332 


Er 


164 


-6.154 


-8.686 


-11.562 


-13.157 


-5.924 


-7.176 


-9.938 


-11.256 




166 


-6.549 


-7.520 


-11.052 


-12.621 


-6.019 


-6.213 


-9.701 


-10.603 




168 


-4.036 


-5.569 


-10.784 


-12.043 


-5.695 


-5.555 


-9.492 


-9.904 




170 


-5.741 


-4.881 


-10.998 


-11.58 


-4.968 


-5.157 


-9.289 


-9.464 


Yb 


164 


-8.651 


-9.462 


-12.069 


-11.256 


-6.388 


-9.336 


-10.435 


-11.708 




166 


-7.919 


-9.114 


-11.527 


-10.603 


-6.459 


-8.223 


-10.427 


-11.350 




168 


-6.702 


-7.947 


-10.946 


-9.904 


-6.024 


-6.715 


-10.12 


-10.771 




170 


-4.353 


-6.061 


-10.603 


-9.464 


-5.274 


-5.105 


-9.795 


-10.051 
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TABLE III. Moments of inertia 2J^^^ in units of MeV~^ for typical well deformed nuclei in 
the rare earth region. The experimental values are extracted from the energies of the first excited 
2"*" states given in Ref. |^^. For other details of table see caption of Table 1. 







without project 


ion 


with projection 


expt. 




A 


RHB 


Gogny 


ixi±D 


Gogny 




Gd 


154 


78.19 


64.79 


/in o/^ 

49.26 


AO C^C 

48.00 


48.75 




156 


88.79 


78.74 


62.15 


67.96 


67.44 




158 


86.38 


79.05 


66.42 


70.90 


75.46 




160 


87.30 


81.56 


r'O on 

68.80 


66.50 


79.72 


Dy 


156 


63.22 


57.52 


4:6.62 


47.79 


43.52 




158 


80.45 


72.81 


56.17 


65.66 


60.64 




160 


93.76 


76.14 


63.71 


67.38 


69.13 




162 


98.22 


77.56 


0/ .z4 


DO.UU 


74.38 




164 


99.92 


82.41 


69.80 


70.69 


81.75 


Er 


164 


88.49 


72.50 


64 7S 


64 62 


65.65 




166 


83.34 


76.40 


68.20 


69.40 


74.46 




168 


93.74 


82.19 


69.04 


69.23 


75.19 




170 


82.30 


81.08 


67.39 


71.55 


76.25 


Yb 


164 


61.66 


60.72 


49.88 


57.10 


48.66 




166 


70.54 


67.72 


61.24 


62.44 


58.61 




168 


82.96 


75.00 


67.37 


68.91 


68.39 




170 


93.12 


80.53 


68.61 


71.69 


71.21 
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Gogny 



. ' 1 ' 1 ' . 

- ^ - 

rf Dy 1 

■ . 1 . 1 . ■ 


. ' 1 ' 1 . 

" ;<:^^©«=£ " 

- • 

■ . 1 . 1 ■ 


. ' 1 ' 1 . 

\ Yb ^ 
■ . 1 . 1 ■ 




1 ' 1 ' 1 

- "^'^^^^^^ " 
■ 1 . 1 . 1 ■ 


1 ' 1 ' 
■ 1 . 1 . " 


1 ' 1 ' 
■ 1 . 1 . ■ 



156 160 164 164 168 164 168 172 

Mass number A 



CRHB 







Gogny 



> -10 
CD 



03 

CO 

"D) 
Q) 

CD 

D) 



-5 - 



CT5 

CL -10 



152 156 160 156 



160 164 164 168 

Mass number A 



. ' 1 ' 1 . 
neutron 


. 1 ' 1 ' 1 . 


. 1 ' 1 ' . 

° CTy^ 

^^^^^^^^^^^^^^^^^^^^ 


" — ^ ^ " 

L Gd J 
■ . 1 . 1 " 


L x""by j 

✓ 

" 1 . 1 . 1 ■ 


" 1 . 1 . ■ 


. ' 1 ' 1 - 


. 1 ' 1 ' 1 . 

' / 


1 1 1 1 1 1 1 1 1 
1 1 1 1 1 1 1 1 1 


^^^^^^^^^^^^^^^^^^^^^^ * 

proton 
■ . 1 . 1 ■ 


■ 1 . 1 . 1 ■ 


■ 1 . 1 . ■ 





1 



1 



164 168 



